Multiple scattering of photons by atomic hyperfine multiplets 
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Mesoscopic interference effects in multiple scattering of photons depend crucially on the internal 
structure of the scatterers. In the present article, we develop the analytical theory of multiple photon 
scattering by cold atoms with arbitrary internal hyperfine multiplets. For a specific application, we 
calculate the enhancement factor of elastic coherent backscattering as a function of detuning from an 
entire hyperfine multiplet of neighboring resonances that cannot be considered isolated. Our theory 
permits to understand why atoms behave differently from classical Rayleigh point-dipole scatterers, 
and how the classical description is recovered for larger but still microscopic objects like molecules 
or clusters. 
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PACS numbers: 42.25.Dd, 32.80.-t, 03.65.Nk 

I. INTRODUCTION 

Atomic physics and quantum optics currently discover 
the fascinating field of mesoscopic physics Q . Mesoscopic 
phenomena are due to interference effects that survive 
a disorder average in phase-coherent samples. Interfer- 
ence effects can be observed using both genuine quan- 
tum matter waves such as electrons or ultra-cold atoms 
and classical waves such as acoustic or electromagnetic 
waves 2J. Arguably, the most dramatic effect is the 
breakdown of diffusive transport due to strong localiza- 
tion, a phenomenon invoked by Anderson in the context 
of the metal-insulator transition 3]. Strong localisation 
has been observed unambiguously for electrons Q and 
microwaves 5] in quasi-one-dimensional systems. Yet 
even in dilute samples far from the dense regime where 
strong localisation could be expected, interference effects 
can be measured. In optics, a rather robust interference 
phenomenon is the coherent backscattering (CBS) effect. 
Here, the constructive interference between waves coun- 
terpropagating along a given multiple scattering path 
enhances the average diffuse intensity reflected from an 
ensemble of random samples in a narrow angular range 
around the backscattering direction 0]. Under optimal 
experimental conditions allowing to apply the reciprocity 
theorem 0, this two- wave interference enhances the in- 
tensity exactly by a factor of 2 @. 

Quantum optical systems involving multiple scattering 
of photons by atoms are well adapted to study general 
concepts of quantum transport. Indeed, atoms are very 
efficient point scatterers for light because the scattering 
cross section close to an internal resonance is huge com- 
pared to the actual atomic size. Laser-cooling techniques 
permit one to prepare low temperature clouds where a 
negligible Doppler broadening of atomic transitions pre- 
serves the phase coherence of the propagating photons. 
For low enough laser intensity, photons are scattered 
completely elastically from closed atomic dipole transi- 
tions such that inelastic scattering and absorption are 



absent (see, however, || for saturation effects at higher 
field intensities). Furthermore, by injecting and detecting 
well-defined photon polarization states, one is capable of 
probing the internal spin degrees of freedom of the atomic 
scatterers. 

The experimental observation of CBS of laser light by 
cold atoms 0, |n| has revealed the crucial impact of 
the internal atomic structure onto coherent photon trans- 
port. Atomic dipole transitions with a nondegenerate 
ground-state scatter photons like isotropic point-dipoles 
(also known as Rayleigh scatterers) 01 • I n general, how- 
ever, photons probe dipole transitions with rather high 
Zeeman degeneracy due to large total angular momen- 
tum. It was shown theoretically that an average over 
all possible angular momentum orientations then leads 
to an antisymmetric component of the atomic scattering 
tensor which reduces the interference contrast of coher- 
ent backscattering considerably 0, Q] . Indeed, in the 
absence of an external magnetic field, photons are scat- 
tered elastically by freely orientable atomic dipoles just 
like electrons propagating in a sample of magnetic impu- 
rities (spin- flip scattering), and in this case interference 
corrections to the diffusive transport picture are cut off 
very efficiently [l5j . 

Previous analytical studies of the CBS double scatter- 
ing contribution from an infinite atomic medium |14| as 
well as the case of multiple scattering inside an infinite 
medium 0] treated an arbitrary single, isolated, degen- 
erate atomic dipole transition. This description is a pri- 
ori accurate if the hyperfine optical resonance under con- 
sideration is sufficiently far away from neighboring transi- 
tions and if the laser detuning is sufficiently small. How- 
ever, Kupriyanov and co-workers |lllll6| showed that the 
CBS enhancement factor displays a slight asymmetry as 
function of the probe frequency around the resonance 
that can only be accounted for by including the other 
optical hyperfine transitions. In the framework of their 
purely numerical calculation, they did not give a qualita- 
tive explanation for the fact that the largest asymmetry 
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is observed in the h\\h polarization channel of preserved 
helicity and that the higher enhancement is found to- 
wards the blue-detuned side of the transition. 

Extending our previous results, we develop in the 
present contribution an analytical theory for multiple 
photon scattering from atoms with a uniform statisti- 
cal distribution over the magnetic quantum numbers in- 
side each hyperfine level of the ground state, leaving 
the possibility of an arbitrary population distribution of 
the various hyperfine levels (for example a thermal one). 
We rely on the decomposition of the scattered intensity 
into irreducible components with respect to the rotation 
group. This method has been developed some time ago in 
atomic physics 0, , but only today, with the advent 
of laser cooling techniques, we can study subtle interfer- 
ence effects for the multiple coherent scattering of quasi- 
monochromatic radiation. As an application, we calcu- 
late the CBS enhancement factor as a function of probe 
frequency for scattering from an infinite half-space. This 
choice is not meant to describe a realistic experiment but 
minimizes the purely geometrical effects of a finite optical 
thickness that varies with the scattering mean-free path. 
Our approach explains sign and magnitude of the asym- 
metry observed in Furthermore, our theory permits 
to describe precisely how the effective degeneracy can be 
reduced by tuning far from hyperfine-structure or fine- 
structure multiplets. This allows us to understand how 
more complicated objects such as molecules or clusters, 
that involve more and more transitions, eventually scat- 
ter light as classical objects with an optimal interference 
contrast. Therefore, we are able to bridge the gap to- 
wards effective theories for mesoscopic photon physics in 
systems with classical scatterers [l9|. 

In Sec. [H] we develop the general theory of multiple 
photon scattering by hyperfine multiplets in the indepen- 
dent scattering approximation. In Sec. lIIII we discuss the 
transition to classical scattering properties. Sec. lIVI con- 
tains the limiting cases of purely clastic scattering and 
the calculation of CBS from atoms with a hyperfine mul- 
tiplet of overlapping resonances. With Sec. [V] we con- 
clude the paper by indicating possible extensions of the 
work. 



II. ANALYTICAL THEORY 
A. Theoretical description 

We describe the photon field for weak laser intensity 
by one-photon Fock states |fee) with wave vector k and 
transverse polarization e. In units where h = c = 1, the 
energy or angular frequency of this state is lo = k. The 
scattering of these one-photon states by atomic dip ole 
transitions will be described following the lines of [141 Hq 
while generalizing to multilevel transitions. 

We consider atoms at rest which are initially prepared 
in an electronic ground-state with angular momentum 
J, but possibly in any sublevel F of the corresponding 
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FIG. 1: (color online) Typical level structure of an atomic 
dipole transition in the L5'-coupling scheme (energy splittings 
not drawn to scale). First the electronic spin S (S = 1/2 in 
this example) is coupled to the orbital angular momentum 
L (here L — in the ground state and L = 1 in the ex- 
cited state) to produce the fine structure angular momentum 
J (here J = 1/2 and J' = 3/2; J" = 1/2 not shown). The 
frequency of the fine structure resonance line is uio (dashed 
arrow) . J is then coupled to the nuclear spin I to give the hy- 
perfine angular momentum F . The allowed hyperfine dipole 
transitions between the levels F g and F^ with frequencies u> eg 
are marked with arrows. 



hyperfine multiplet. The probe light frequency is near- 
resonant with an electric dipole transition connecting the 
electronic ground-state J to an electronic excited state 
,/' = J — 1, J, J + 1. We assume here that the atom 
fine structure is well described by an LS'-coupling scheme 
(Fig. ^) , where the total electronic angular momentum 
J = L + S couples to the nuclear spin I, producing a 
hyperfine multiplet F = J + I with F = \I — J\, .... I + J . 
The energy splitting of each hyperfine multiplet then is 
much smaller than the fine structure splitting. This is 
the typical situation in alkali atoms like Rubidium: for 
the D2 line (5S 1/2 5P 3/2 ) of the Rb 85 isotope with 
nuclear spin / = 5/2, the ground-state J = 1/2 is split 
into two levels with F = 2 and F = 3. Similarly, the 
excited state J' = 3/2 is split into F' = 1,2,3,4. In 
this scheme, the effect of different far off-resonance fine- 
structure transitions J — > J" with the same hyperfine 
quantum numbers F — > F' is neglected. Other coupling 
cases could be easily treated along the lines of reasoning 
described in the present paper. In general, several opti- 
cal transitions F — > F' = F — 1, F, F + 1 between the 
ground-state and the excited state hyperfine multiplets 
are allowed. We will assume that the J — > J' transition 
is closed so that, after having scattered an incoming pho- 
ton, the atom returns to the same ground-state hyperfine 
multiplet but possibly in a different hyperfine level (this 
is the case for the Rb D2 line). 

The Hamiltonian of this internal atomic structure with 
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ground and excited levels g and e is 



(1) 



where w 9 is the energy of the atomic level with angular 
momentum F g . P g is the projector onto that level, 



P 9 = E \ F 9 m 9)(F 9 



(2) 



with a similar expression for the projector P' e onto excited 
levels with angular momentum F' e . Throughout the pa- 
per, the following convention holds: primed symbols like 
F' e refer to the excited state multiplet whereas unprimcd 
symbols like F g refer to the ground-state multiplet. 



B. Scattering amplitude 



natural linewidth T. Note however that this approxima- 
tion can easily be relaxed: in the formulas below describ- 
ing interference effects between various hyperfine com- 
ponents, such as Eq. i|f 9[l . one would have to use the 
frequency-dependent T eg (u;). Since the dipole operator 
D acts only on the electronic angular momentum ,7, and 
not on the nuclear spin J, the reduced matrix element 
can then be reduced even further [20|: 



((J'I)F'\\D\\(JI)F) =(-) J ' +I+F+ W(2F + l)(2F< + 1) 

J' F' I 
F J 1 

(6) 



(J'\\D\\J) 



The 6j-symbol describes the recoupling of the 4 angular 
momenta I, J, F, and 1 for the vector operator D. The 
sum over the final hyperfine levels Ff in |5| then can be 
done using a 6j-symbol normalization rule such that one 
finds a uniform width for all excited hyperfine levels, 



While scattering a single photon \ke) i— > \ke), the atom 
undergoes a transition F g — > Ff from an initial to a final 
hyperfine level. The corresponding scattering amplitude 
T g f(ke i— > fee) = Y^ij ^t^ij { L0 ) e j is a matrix element of 
the scattering operator T 9 ^ (uj) acting on internal atomic 
states and polarization vectors. According to fundamen- 
tal rules of quantum theory, the amplitudes for all indis- 
tinguishable scattering processes via intermediate excited 
states e have to be added coherently: 



(3) 



The Cartesian components of the scattering operator for 
the partial scattering process F g — > F' e — ► Ft are [Tj| 
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2eoV uj + uj g — uj e + iT eg (uj)/2 



(4) 



Here D is the electronic dipole operator with components 
Di. Because of energy conservation the final photon fre- 
quency is ujf g = u) + ujg — u>f. The scattering is elastic 
when LUf g — lu, i.e., when the initial and final state have 
the same energy u> g = u>f. It is inelastic in all other 
instances. 

The decay rate T eg (uj) of each excited level is the sum 
of spontaneous decay rates to all accessible final ground- 
states Ff. 



T eg (u) 
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For atoms used in cold atom experiments such as alkali 
or alkaline earth atoms, the hyperfine splitting is much 
smaller than uj such that the frequency difference between 
the various hyperfine transitions can be safely ignored for 
the evaluation of T eg (uj), which is then a constant, the 



r = 



37Te 



(7) 



Here, we introduce the reduced fine structure matrix 
element d = (J'\\D\\J)\/ y/2J' + 1. For further use, we 
define the dimensionless dipole operator d eg = P' e DP g /d 
between the levels F e and F g with detuning 

5 eg — uj + uj g - u e . (8) 

Finally, the scattering operator takes the form 



r/2 



2vrp(w) S eg +iT/2 



(df e )i(d eg )j 



(9) 



The prefactor contains the free photon spectral density 
p{uj) — Vui 2 /2ir 2 for a given polarization in a quantiza- 
tion volume V. 



C. Averaging over atomic degrees of freedom 

Multiple scattering in the dilute regime can be depicted 
as a succession of scattering events by single atoms con- 
nected by propagation in an effective medium. An effec- 
tive photon transport theory is obtained by a configura- 
tion average over atomic degrees of freedom. The atoms 
are assumed to be initially prepared independently in the 
ground-state levels F g with probabilities p g . These prob- 
abilities can represent an equilibrium distribution at a 
certain temperature or situations where hyperfine pump- 
ing has been achieved. Inside each level, we assume a 
uniform statistical distribution over magnetic quantum 
numbers m g . Accordingly, the one- atom density matrix 
of internal degrees of freedom reads 



Pat 



E 



P9 4 0) 



2F g + l 



(10) 
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The assumption of a scalar pf^ is reasonable for opti- 
cally thick samples, where the isotropization by multiply 
scattered photons dominates possible optical pumping ef- 
fects due to the incident laser. Needless to say, this as- 
sumption greatly simplifies the calculation. In essence, 
an entirely analytical description is only manageable be- 
cause the average over an isotropic distribution like Hl()|) 
restores rotational invariance. 



D. Scattering mean free path and total cross 
section 

The wave vector k of a photon with frequency u in the 
dilute atomic medium is determined by the dispersion re- 
lation k(u>) = Lon r (u>) where n T (ui) is the refractive index. 
In scattering media, the refractive index has an imagi- 
nary part that reflects the fact that scattering depletes 
the propagating average photon field. This defines the 
scattering mean free path Imfc(cj) = l/2£(u>). Techni- 
cally, the scattering mean free path is calculated via the 
photon self-energy, which in turn is proportional to the 
scattering o pera tor averaged over the atomic density 
matrix i|10|) 14, 15]. Carrying out the isotropic sum over 
m g , one finds 



— (2j+i) ^ Pg i^YT 



(11) 



Here, n is the number density of atoms. The sum over 
the ground-state levels g is due to the simple atomic dis- 
tribution IjlUfl. The sum over the excited levels e is due 
to the linear superposition of scattering amplitudes, and 
the Lorentzian frequency dependence comes from taking 
the imaginary part of the resonant denominators in J2J . 
The atomic level structure enters via the coefficients C eff , 
a short-hand notation for the 6j-symbols 



J' F' e I 
F a J 1 



(12) 



that stem from the decomposition 10. Thanks to the 
6j-symbol normalization, they obey the following sum 
rules E3: 
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^2(2*1 + 1) Ct; 
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E(^ + i)^ 



2J + 1 
1 



(13a) 
(13b) 



Energy conservation dictates that extinction of the for- 
ward propagating mode is related to the total scattering 
cross-section cr(u>) by the relation ^(w)^ 1 = ncr(cj). This 
relation, the optical theorem in a multiple scattering dis- 
guise, will be proved in the following subsection. We can 
use it here to write the total scattering cross-section as 



l) C e S 



4<5 e yr 2 



(14) 



This total cross section is a weighted sum of Lorentzians 
of uniform width for all possible transitions (g — > e — > 
/). Note that this result pertains to scattering from a 
single hyperfine multiplct in both the ground and the 
excited states, i.e., a single fine-structure transition J — > 
J' . Moreover, we consider a single dipole transition with 
a single principal quantum number in the excited state 
such that no interference effects between different fine- 
structure transitions or different dipole transitions with 
the same angular momentum (as considered in \22\) can 
appear in the total cross section. 

Since the total cross section measures the total deple- 
tion of the initial photon state, it includes both elastic 
scattering events where the final atomic level / is iden- 
tical to the initial level g, but also inelastic scattering 
where the internal energy of the atom is changed. A 
separate evaluation of elastic and inelastic scattering is 
possible by considering the differential cross section. 



E. Differential cross section 

In the weak scattering approximation, the building 
block of multiple scattering is the average differential 
cross section for scattering polarized photons \ke) i— ► |fee) 
while the atom undergoes a transitions (g — ► /). The dif- 
ferential cross-section is obtained by averaging the square 
of the scattering operator © over the atomic density 
matrix (|10H . giving the well-known Kramers-Heisenberg 
formula [23j . When all scattering processes without final 
frequency analysis are considered, it can be written in 
the form 



da 



i<j{<jj) 
8tt 



x(e, e , e, e 



(15) 



Here, the total cross section cr(uS) has been factorized 
such that the dimensionless atomic vertex funtion I con- 
tains all angular information about light polarization and 
internal atomic structure. This vertex function has been 
calculated in 0] f° r the case of a single isolated reso- 
nance J — * J' using the techniques of irreducible tensor 
operators. Since we still assume a scalar density ma- 
trix (|ll)p. the vertex function must still be the sum of all 
possible scalar products of its vector arguments, 

X(e, e* , e, e*; u) — W\ \e* ■ e| 2 + w 2 |e • e| 2 + w 3 . (16) 

As in 0] i the weights are given by the combinations 



Wi 



so - s 2 



w 2 



»2 



W 3 



s 2 + si 



(17) 



of coefficients Sk that describe scalar (A" = 0), antisym- 
metric (K = 1), and symmetric (K = 2) scattering. For 
the present case of multiple resonances, these coefficients 
are now frequency dependent: 



s K (u)) = J2p9J2 s k^ 

9 f 



(18) 
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where each transition g — > / is described by 



F. Elastic ds. inelastic scattering 



K 



()7T 



■(2J' + 1) 2 (2 J F / + 1) 



e/ 



i - 2W eg /r 



k 2 a{uj) 

The structure of these expressions reflects well-known 
rules of quantum theory: the total transition probability 
is the weighted incoherent sum over initial states g with 
occupation probability p g , and the incoherent sum over 
distinguishable final states /, but the square of the coher- 
ent sum of all indistinguishable amplitudes, here related 
to the intermediate excited levels e. The atomic level 
structure enters via the coefficients C eg and C e f defined 
in l|12|l . The essential information about the irreducible 
tensor modes K = 0,1,2 of the scattering operator is 
contained in the factors 



*4 9/ = (-l)^(2F e ' + lh F p pi 



1 1 K 



(20) 



The coefficients sf/(w), together with the scattering 
mean free path (|lll) . are the basic ingredients for com- 
puting multiple scattering quantities. The derivation of 
expressions <|18|) - (I20I1 constitutes the main achievement of 
the present work. From this point on, we will essentially 
explore its consequences. 

The total cross section (|14|l can be calculated from 
(|f 5f) by a sum over all final photon polarization vectors 
e and an angular integration over the scattered photon's 
direction fc/|fc|. This operation on the vertex function 
(|Td|) yields (8tt/3) x ( Wl + w 2 + 3w 3 ). The weights (H7J) 
have been defined such that they obey the sum rule |14| 

wi(w) + w 2 {uj) + 3w 3 (uj) = 1. (2f) 

which is equivalent to 

(2K+l)a K (w) = 3. (22) 

This last relation can be deduced from l|18|l by virtue of 
6j-symbol orthogonality |2l| and the sum rule (|13b|l . In a 
more general setting, this relation is shown to be a trace- 
conservation property of the intensity scattering vertex 
[24|. With this, we indeed recover the total cross section 
(|I4|I and thus prove the optical theorem £(ui) = l/na(w) 
that was used above to derive the total scattering cross 
section directly from the mean free path. 

Note that the interference between scattering ampli- 
tudes via different excited states e, present in (|I8|I . disap- 
pears in the total cross section <|14|) . This is a result of the 
complete statistical average over the degenerate ground- 
state. On the elementary level of Clebsch-Gordan coef- 
ficients, the somewhat abstract 6j-symbol orthogonality 
appears as a complete cancellation of interference terms 
from equiprobable ground-states for large detuning. 



Expression l|18|) of the differential cross-section coeffi- 
cients permits to distinguish between elastic scattering 
events with f — g and inelastic scattering events with 
/ 7^ g such that we are able to write 



da \ - x - /dcr\ 



(23) 



where each transition g — ► / is described by its coeffi- 
cients s 9 J (lu) defined in l|19|) . By integrating over final 
polarization and scattering directions, we can therefore 
write 

a(uj) = CTci(cj) + cr incl (cj) = ^Pg^ (J 9 E f^ UJ ) ( 24 ) 

9 f,e 

with the total cross section for each elementary transition 
9 -*■ e -> f, 

^ / M = |(2J' + I) 2 (2F / + l)^f^^. (25) 

By separating elastic from inelastic contributions, we 
therefore find cr(w) = a c \(uj) + cri ne i(u;) with 



(26) 



O-inclM = J2 P 9 CT 9 e /H- ( 27 ) 

9 f^9,e 

In principle, a frequency analysis of the scattered photons 
makes it possible to measure elastic and inelastic contri- 
butions independently. Note, however, that all elastic 
components from the various initial states g are at the 
same frequency and cannot be distinguished by a fre- 
quency analysis. 

Elastic scattering is of course desirable for coherent 
multiple scattering since the scattered photons stay on 
resonance. Therefore, previous experiments of coherent 
backscattering from degenerate atomic dipole transitions 
have been performed on various isolated hyperfine tran- 
sitions F — > F' that are closed (see (2f| for a compre- 
hensive list). Or rather, these transitions are very nearly 
closed: a small rate of inelastic scattering persists be- 
cause after excitation to an off-resonant hyperfine level, 
inelastic transitions are possible, for example the process 
F g = 3 — ► Fg = 3 — * Ff = 2 is the case of the much- 
studied transition F = 3 1 — > F' — 4 of Rb 85 . Detuning 
away from the closed resonance towards the other reso- 
nances of course increases the inelastic scattering rate. 
A completely consistent theory of scattering from hyper- 
fine multiplets therefore needs to take into account in- 
elastic scattering, without which predictions about CBS 
enhancement factors close to open transitions are ques- 
tionable Inelastic scattering can be incorporated 
in Monte Carlo simulations of photon trajectories but is 
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beyond the scope of the present article, devoted to an- 
alytical results. In Sec. II VI we will therefore specialize 
to the interesting case where hyperfine multiplets are si- 
multaneously excited from a unique ground state which 
assures purely elastic scattering. 



III. TRANSITION TO CLASSICAL 
SCATTERING PROPERTIES 

The complicated radiation pattern described by 
Eq. I jltjjl is generally not the one of a classical point- 
dipole scatterer. But for an isotropic atomic transition 
F = — * F' = 1, there is only a single nondegenerate 
ground state and a single threefold degenerate excited 
state. In such a case, the weights coefficients in (|16H are 
simply wi — 1,102 = if 3 = 0, and the atom scatters light 
like a classical point-dipole or Rayleigh scatterer with 



da 



8tt 



|e*- e | 2 , a(u) = 



()7T 



1 



fc 2 i + 4<5 2 /r 2 



(28) 



This is equivalent to the quasiclassical model of an elas- 
tically bound electron. In this case interference effects 
in multiple scattering are fully preserved and the CBS 
enhancement factor achieves its maximal value 2 in the 
helicity-preserving channel [12| . 

In optics, scattering of light — including interference ef- 
fects like CBS — by small particles (smaller than the wave- 
length of light) is often successfully described by mod- 
elling the small particles as classical point-dipole scatter- 
ers. However, such classical particles are very unlikely to 
be in a pure F = state, and we know that atoms in 
F > states show ver y lo w CBS interference in all po- 
larization channels (lot fl3| . Therefore, the validity of the 
classical point-scatterer model for complex objects with 
possibly many internal resonances is rather surprising. 
In this section, we therefore show how "classical" coher- 
ent backscattering properties are recovered for complex 
quantum objects (like atoms, molecules, or clusters) us- 
ing a time-scale argument together with the analytical 
formulation developed above. In other words, we close 
the gap between the microscopic theory for atoms and 
effective models for classical scatterers. 



A. Time scales 

The fact that hyperfine and fine structures affect the 
angular distribution of the scattered light has been recog- 
nized and understood for quite some time [2||. In order 
to grasp how coherent light scattering depends on the 
hyperfine structure (or fine structure) multiplets and the 
detuning from resonance, it is useful to consider the var- 
ious time scales coming into play. 

Since F is a good quantum number, the semiclassi- 
cal picture of hyperfine structure is that both the nu- 
clear spin J and the electronic angular momentum J 



precess around their constant sum F = I + J . It fol- 
lows that the atomic electric dipole precesses after being 
excited by the incoming laser field. If the laser source is 
very monochromatic — lincwidth smaller than the classi- 
cal precession frequency, i.e., temporal coherence longer 
than the precession period — the radiated field is built by 
a coherent superposition of radiating dipoles with various 
spatial orientations. The net result of this coherent su- 
perposition, calculated using quantum theory of angular 
momenta, is the unusual radiation pattern There 
are, however, two complications: first, the atomic dipole 
decays because of spontaneous emission over the time 
r _1 ; second, if the excitation is not exactly resonant, it 
also oscillates at a frequency equal to the detuning 5 from 
resonance. 

These two effects are properly understood by consider- 
ing the scattering of a quasimonochromatic wave packet 
formed by superposition of different neighboring frequen- 
cies around a central frequency u). Such a wave packet 
is not scattered instantaneously, but after some delay 
tw(w), known as the Wigner time delay. It is given by the 
derivative of the scattering phase shift arctan(2<5/r) — see 
@ — with respect to frequency: 



*w(w) 



2 

f r 



i 



4<5 2 /r 2 



(29) 



The Wigner time delay is the time scale after which the 
atomic dipole induced by the incoming electromagnetic 
field ceases to radiate coherently. It is maximum at reso- 
nance where it is twice the lifetime of the atomic excited 
state, and decays towards zero away from resonance. 

If the Wigner time delay iw(w) is longer than the pe- 
riod of hyperfine precession Ta = 27r/A, where A is the 
typical hyperfine splitting, then the radiating dipole pre- 
cesses during its coherence time, thus giving rise to a 
specific radiation pattern. This is the case at resonance 
if the linewidth T is smaller than the hyperfine splitting 
A, a situation usually encountered for an isolated hyper- 
fine component. By contrast, if the Wigner time delay 
t-w(u>) is shorter than the period of hyperfine precession 
Ta, then one expects the radiating dipole to be spatially 
frozen during its coherence time. Consequently, the ra- 
diation pattern should turn into the one associated with 
a classical dipole. 

There are two obvious ways of reaching the classical 
scattering situation £w( w ) T\. On the one hand, if 
the spontaneous decay rate is larger than the hyperfine 
splitting, r > A, the various components of the hyper- 
fine multiplet are not resolved. The degree of freedom 
associated with the nuclear spin and its interaction with 
the electron are completely quenched on the radiative 
decay time scale and can be forgotten, whatever the de- 
tuning from resonance. On the other hand, when the 
resonances are well separated, T <C A, one can use large 
detuning to bring the Wigner time delay below the hy- 
perfine precession period. By detuning far away from 
the entire multiplet of hyperfine resonances, the excita- 
tion then no longer probes the nuclear spin, and therefore 
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scattering is only sensitive to the fine structure. 



B. Analytical derivation 

The full analytical theory developed in Sec. ITU permits 
one to follow precisely how the different resonance con- 
tributions combine to yield classical scattering character- 
istics for large detuning. When the Wigner time delay is 
shorter than the hyperfine period, all resonant denomi- 
nators (5 e g +i r/2) can be taken equal to a common value 
(6 + iT/2), where 5 = to — loq is now the detuning from 
the fine-structure resonance line J — > J' (cf. Fig. [JJ. In 
expression (|14|l . the resonant denominator can thus be 
factorized from the sum over excited states e. The sum 
rule l|13a() together with the normalization J2 g Pg = 1 can 
then be used to obtain the effective fine-structure cross 
section 



where 



Mjj, = 



6tt M 



i.v 



i + 4<P/r 2 



2 J' + 1 
3(2 J + 1) 



(30) 



(31) 



is the ratio of multiplicities normalized to Moi = 1. The 
cross section is now independent of the hyperfine struc- 
ture, i.e., no longer depends on the nuclear spin / and the 
population p g of the various hyperfine levels. Similar ar- 
guments apply to the scattering of broadband radiation 
by thermal atoms (see especially Sec. 3.4.2. of |l7j)- 

The radiation pattern of the scattered photons is 
subject to the same transformation: the frequency- 
dependent denominator can be factorized from the sum 
over e in l|18Jl . For large enough detuning, the atomic 
medium makes no difference between photons scattered 
clastically or inelastically, such that all contributions of 
final states / must be added. The resulting sum involves 
the product of three 6j coefficients which can be com- 
puted using the Biedenharn-Elliott sum rule |2lJ and ev- 
erything boils down to the coefficients calculated in 



s K = 3(2J' + 1) 



1 1 K 
J J J' 



(32) 



The differential cross section then is still given by l|15|l. 
its only frequency dependence being the single Lorentzian 
resonance of the total cross section ()30|l . 

The very same arguments apply to the electron spin 
responsible for the atomic fine structure. Indeed, the 
coupling of the electronic orbital angular momentum L 
with the optically inactive electronic spin S produces the 
total electronic angular momentum J = L + S . Again, 
when the spontaneous decay rate Y is larger than the fine- 
structure splitting (meaning that the various components 
of the fine-structure multiplet are not resolved) , or when 
the laser frequency is far detuned from the fine-structure 



multiplet such that one has to sum over all possible ex- 
cited levels J', one recovers the case of a L — ► L' tran- 
sition, where only the orbital properties (directly related 
to the charge density response to the laser excitation) 
of the electrons play a role. The formulas arc simply ob- 
tained from eqs. l]3Up. I|31|) . and i|32[l through replacement 
of (J, J') by (£,£'). 

If the ground state of an atom is an S state with L = 0, 
which can be optically excited only to a P state with 
L' = 1, we find Mll< = 1 and the sk coefficients (|32f) 
simply reduce to Sq = 3, s\ = S2 = 0. We are back to 
the situation i|28|) where the atom radiates exactly like a 
classical point-dipole scatterer since only a single atomic 
transition — > 1 is optically active. 



C. Molecules and complex objects 

The previous analysis can be extended to objects 
slightly more complex than atoms, but whose energy 
spectrum and eigenstates can still be calculated. Let us 
for example consider a diatomic molecule with rovibra- 
tional structure. If the molecular linewidths are suffi- 
ciently small for the rotational structure to be resolved, 
and if the incoming light is sufficiently monochromatic for 
a single rotational line to be resonant, then the molecule 
will scatter light with the specific radiation pattern of the 
resonant J — > J' transition. This is because the Wigner 
time delay is longer than the rotational period of the 
molecule. In contrast, if several J — » J' = J, J ± 1 ro- 
tational transitions have to be taken into account coher- 
ently, a sum over excited states has to be performed, very 
similar to the one in @, with the only difference that 
the dipole matrix elements depend now on the molecu- 
lar quantum numbers, instead of J and F. These dipole 
elements are given in [27| for various possible couplings. 
A full discussion is beyond the scope of this paper, but 
the net result is as expected: the sum simplifies thanks 
to sum rules over 3j and 6j symbols and the radiation 
pattern depends only on the electronic dipole transition 
considered. In the most common case where it is a S — > II 
transition (the molecular equivalent of an atomic S — ► P 
transition), one recovers again as a net result the fact 
that the molecule scatters light like a Rayleigh scatterer. 

However, when several electronic states have to be 
taken into account — either in atoms or in molecules — 
interference between several transitions is in principle 
still present. When one considers larger objects, such 
as polyatomic molecules for example, the time scales are 
considerably affected. The total moment of inertia in- 
creases, making the rotation slower and the correspond- 
ing classical rotation time longer. At the same time, ex- 
citations of the object decay faster, making the Wigner 
time delay shorter. Rotational degrees of freedom are 
thus quenched. Whenever the internal evolution time 
is longer than the Wigner time delay, one expects the 
internal structure to be frozen and irrelevant for light 
scattering properties. Generically, the shortest internal 
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time scale which remains shorter than the Wigner time 
delay, will be associated with the electronic excitation, 
which is the degree of freedom excited in the optical fre- 
quency range. An object with more than few atoms will 
thus scatter "classically" , i.e., following the average elec- 
tronic polarizability [28| . For an object much smaller 
than the light wavelength (such as a cluster with, say, 
1000 atoms), the multipolar expansion of the electronic 
polarizability will be dominated by the dipolar term, pro- 
ducing the dipole radiation pattern (J2BJ. There may of 
course be cases where the dipole contribution vanishes 
and higher orders dominate the radiation pattern, but 
such situations are certainly not generic. Finally, for 
objects larger than the light wavelength, the shape of 
the object is resolved and determines the light scattering 
properties. When the scatterers are anisotropic and pos- 
sibly oriented, for example nematic liquid crystals in ex- 
ternal fields, the scattering properties remain anisotropic 
and are described by a classical anisotropic polarizability 
tensor 291. 



IV. ELASTIC SCATTERING 



J' 



-F' 

: F' 

mi 

- F 

J- m, 

- F 



J' = l 



J = 



F'=I+1 
-F'=I 
•F'=I-1 



-F = J 



(a) 



(b) 



(c) 



FIG. 2: (color online) Level scheme for limiting cases of purely 
elastic light scattering, requiring a unique ground-state F g . 
Effectively excited levels are enclosed by a dotted circle, (a) 
Zero nuclear spin: only the fine structure with a unique ex- 
cited state is involved, (b) Resonant excitation of a closed 
transition F max — > F max + 1 or F min — » F min - 1, also with a 
unique excited level, (c) J = 0: a unique ground level with 
F — I is coherently coupled to several excited levels. Cases 
(a) and (b) reduce to the known results of 14]. Case (c) 
with a structured excited multiplet is treated analytically in 
section IIVI 



The general expressions of Sec. ^simplify considerably 
in the case of a unique ground state since only elastic 
scattering can occur. This is the case when there is no 
hyperfine splitting (because either / = or J = 0), or 
when a closed dipole transition far from other resonances 
is probed (cf. Fig. [2J ■ In such cases, the 6j-coefficients 
(|12|l obey certain selection rules and are most easily cal- 
culated using the sum rule eq. I|13b|l that reduces to a 
single term for the only ground level F g : 



C 



1 



(2J' + 1) (2F g + l)' 



(33) 



an expression which is notably independent of the excited 
state(s) e. In the following, we distinguish between cases 
where a single excited state is relevant (and we are back 
to the previous results and the new situation where 
several excited states have to be considered. 



A. Unique excited state 

In the two situations shown in Fig.^a) and (b), there 
is only a single effective two-level system to be considered. 
In these two cases, sums over ground and excited states 
reduce to a single term and we recover the results of : 

(i) Zero nuclear spin 1 = 0. In this case, there is no 
hyperfine splitting and F g = J, F' e = J'. The coupling 
coefficient J23J) is C 2 eg = 1/(2 J+1)(2J'+1). The net effect 
is the same as for a large detuning from the multiplet 
discussed in Sec. IIIII we recover the total cross section 
(|30|l and the intensity vertex coefficients l|32[l for a J — > J' 
transition. 

(m) On-resonant excitation of a closed transition 



F min = Fmin — 1- The latter case is possible only if 
\I — J\ > 1. Naturally, the hyperfine structure splitting 
A should be large enough such that other resonances can 
indeed be neglected. The total cross section involves only 
a single term and is given by 



, v 6tt M eg 



k 2 1 + 4<s 2 „ /r 2 ' 



with the short-hand notation 



M eg = M FgF , 



2^ + 1 
3(2F S + 1) 



(34) 



(35) 



for the ratio of multiplicities. The differential cross- 
section coefficients (|32() are 



.* = 3(2*2 + 1){^ l pg K F ,) 



(36) 



T — > F' 



1 Or F„ 



A 



These expressions are identical to eqs. (|3TJ|) and 
with substitution of J b y F„ and J' by F' e which was the 
situation anticipated in |14[ . 



B. Multiple excited states 

A more interesting situation occurs when there is a 
unique ground state coupled to a multiplet of excited 
states, e.g., for a vanishing electronic ground-state angu- 
lar momentum J = like in Fig. Now the ground 
level F g — I is unique, and the excited level is split into 
the hyperfine levels F' e = 1—1,1,1+1 (for I > 1) or 
F' e = 1,1+ 1 (for I = 1/2). Using eq. iJSJ, the total 
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cross section (|14fl can be written as 



•4*2/1* 



(37) 



a sum of Lorentzians weighted by the multiplicities H35|) 
of the various excited states. More interesting, and part 
of the central results of this paper, are the frequency- 
dependent intensity coefficients 



6ir(2F g 



k 2 a(uj) 



K 



(38) 



where u e J? 9 is obtained from l|20|) by putting g = f. These 
coefficients permit to describe single and multiple scat- 
tering for arbitrary detuning between and outside the 
resonances. 

Expression (|38|l shows clearly that the coherent super- 
positon of amplitudes © carries through to the aver- 
age intensity: the coefficients sk are squares of interfer- 
ing amplitudes and not the sum of squared amplitudes. 
Therefore, the influence of other resonances may lead to 
subtle phenomena in multiple scattering, sensitive to the 
differential cross section, which are not immediately vis- 
ible in the total scattering cross section l|37l) . Indeed, 
by virtue of a 6j-symbol orthogonality, the total cross 
section is just the "incoherent" sum of the individual 
cross sections and therefore insensitive to these interfer- 
ence effects. In the classical picture, this can be nicely 
understood: the three hyperfine matrix elements actu- 
ally originate from the same optically active transition 
J = — * J' = 1, and the precession of J around the 
nuclear spin J modifies only its spatial repartition, not 
the total scattering rate. 

In elastic coherent backscattering of light by 

atoms has been calculated analytically in the double scat- 
tering approximation and for a semi-infinite scattering 
medium, at fixed detuning from an isolated resonance 
J — > J'. Now, all results of Q can be extended to arbi- 
trary values of light frequen cy. Similarly, using the Monte 
Carlo method described in |3(j], the full CBS cone (with 
arbitrarily large scattering orders) can be computed nu- 
merically. 

As a specific example, we will consider the case of 
the fermionic isotope 87 Sr of strontium where the energy 
splitting between hyperfine components is comparable to 
the widths of the resonances themselves [3l| . Fig. Ofa) 
shows the total scattering cross section for the 1 So — > *Pi 
optical transition (T/2ir = 32 MHz) with F g = I = 9/2 
and Fg = 9/2, 11/2, 7/2. The frequency separation be- 
tween the (9/2 -> 9/2) and the (9/2 -►11/2) resonances 
is 17 MHz = 0.531 T/2n while it is 60 MHz = 1.875 T/2n 
between the (9/2 -> 9/2) and the (9/2 -> 7/2) reso- 
nances. The vertical dotted lines indicate the position of 
each hyperfine resonance which clearly cannot be consid- 
ered as isolated. 




-10 12 3 

detuning S (units of T) 

FIG. 3: (Color online) (a) Total single scattering cross sec- 
tion, eq. 13711 . for Sr 87 as a function of detuning (in units of 



r). The ground-state angular momentum is F g 



9/2. 



There are three accessible excited states with angular momen- 
tum F' e = 9/2, 11/2, 7/2 and frequency separation 0.531 F/2tt 
and 1.875 F/27T, respectively. The corresponding (overlap- 
ping) optical resonances are shown by dashed lines, (b) Fre- 
quency dependence of the single scattering bistatic coefficient 
71 for a semi-infinite medium in the four polarization channels 
ft || ft, ft -L ft, 1 1| I and I -LI. (c) Frequency dependence of the 
double scattering CBS interference contrast f° r the 

same situation. The predictions for isolated resonances [l4ll . 
indicated by the symbols, are clearly off the lines, indicating 
the importance of interference, (d) Frequency dependence 
of the total CBS enhancement factor 1401 1. Monte Carlo cal- 
culation including all orders of scattering for a semi-infinite 
medium. 
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C. Single backscattering from a semi-infinite 
medium 

The amount of light reflected off the sample after N 
scattering events can be quantified by a dimensionless 
parameter called the bistatic coefficient 7jv. The sin- 
gle scattering bistatic coefficient 71 in the backward di- 
rection for a semi-infinite medium with constant spatial 
density exposed to an incoming plane wave at normal in- 
cidence is essentially the atomic polarization vertex 1|16|) 
(see eq. (42) of [l4|, namely 

7i(w) = ^[w 1 (uj)\i*-e\ 2 +w 2 {cu)\i-e\ 2 + w 3 (uj)] . (39) 

In Fig. OJb), 71 is plotted as function of the detuning 
from the transition 9/2 — > 9/2 by using the frequency- 
dependent weights Wi(uj) given by l|17f) in terms of (|38|l . 
The four curves correspond to the usual polarization 
channels: h\\h (parallel helicitics, i.e., opposite circular 
polarizations in the backward direction), /i_L/i (orthog- 
onal hclicities), l\\l (parallel linear polarizations), and 
I _L / (orthogonal linear polarizations) . Since we take a 
semi-infinite medium, every photon entering the medium 
must eventually exit. The atomic internal structure only 
redistributes photons into different polarisation chan- 
nels. The total intensity is independent of frequency, 

j 1 (i\\i) + 11 (i±i) = 7l {h\\h) + ^(h±h) = 3/4. 

Far from resonance, the atom radiates like a point- 
dipole scatterer, in agreement with the general discussion 
of the transition to quasi-classical scattering fsection lTTTll : 
since the radiated polarization is equal to the incoming 
one, the bistatic coefficient vanishes in the I _L I and h\\ h 
channels and takes its maximal value 3/4 in the I \\ I and 
h _L h channels. Close to resonance, the situation is more 
complex. Anisotropic scattering populates all polariza- 
tion channels and interference effects between the various 
hyperfine components are clearly made visible by the po- 
larization analysis. 

D. Double scattering interference contrast 

The multiply scattered intensity (N > 1) contains 
two dominant contributions: first, the so-called ladder 
contribution 7^ without interference between multi- 
ple scattering amplitudes. This corresponds to waves 
co-propagating along scattering paths. Second, the so- 

(C) 

called maximally crossed contribution 7^ which origi- 
nates from counter-propagating waves and incorporates 
the interference responsible for the CBS peak. The CBS 
interference contrast at scattering order N is the ratio 
7]^ ] H^n' For double scattering, N = 2, the formulas 

(B6) and (B24) of [3 directly give 7^ and ] at exact 
backscattering in terms of the weights (u>i, u>2, W3). In 

Fig. OIc), the double scattering CBS contrast 
is plotted as function of detuning. The on-resonant pre- 
dictions for well separated resonances, shown as symbols, 



neglect the interference between the different resonances 
and are clearly not suited for quantitative predictions. 
Outside the multiplet, the double scattering CBS inter- 
ference contrast approaches the maximum value 1 pre- 
dicted for the isotropic dipole or non-degenerate transi- 
tion J — — > J' = 1 in all channels (beyond double 
scattering, the contrast is unity only in the parallel chan- 
nels h || h and 1 1| I), again in agreement with the general 
discussion in Sec. IIIII 

Considering only the resonance with highest frequency, 
here 9/2 — * 7/2, it is evident that the interference con- 
trast is generically larger towards the blue side in order 
to reach the optimal value 1 at large detuning. Further- 
more, the largest asymmetry or slope on resonance is 
found for the h\\h channel since it starts from the low- 
est contrast on resonance. These features, present in the 
figures of 0| but unexplained by the authors, therefore 
find a natural explanation. 



E. CBS enhancement 

The enhancement factor of the CBS cone is 

a = 1 + 1C , (40) 
1L + 7i 

where j L = Ew>2 7jv and 7c = EiV^T^f 5 account 
for all multiple scattering orders. The presently derived 
weights u>i(u>) can be used for evaluating the multiply 
scattered intensity via the propagation eigenvalues of the 
ladder and crossed series, X(lu) and x( w )i as derived in 
[TH- For third order scattering and beyond, exact an- 
alytic calculations become very complicated, and it is 
more convenient to turn to a Monte Carlo approach, as 
described in [H |H H3 . In Fig.^d), we plot the CBS 
enhancement factor for 87 Sr as a function of detuning, 
in the four polarization channels and for a semi-infinite 
medium. The most obvious observation is that the en- 
hancement factors are typically small in the region of 
overlapping atomic transitions, and take larger values at 
large detuning, as expected from the transition to clas- 
sical scattering behavior discussed in Sec. IIIII High or- 
ders of scattering contribute significantly to the ladder 
intensity, but only weakly to the CBS contribution be- 
cause phase coherence is rapidly lost after the average 
over atomic degrees of freedom. The enhancement factor 
thus behaves similarly to the double scattering contrast, 
only amplifying the changes with the detuning. 

A strong reduction of the CBS enhancement factor in 
the vicinity of a particular resonance may have various 
origins. For example, in the h\\h channel, for negative 
detuning, the double scattering contrast is rather high. 
This is because the dominant contribution comes from 
the 9/2 — > 9/2 transition of the J — ► J type with large 
J (see EH). The reduction of the enhancement factor 
here must be attributed to single scattering. For posi- 
tive detuning, the situation is opposite: single scattering 
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is rather low, but also the double scattering contrast is 
poor. 

In the orthogonal polarization channels I _L I and h _L h, 
the multiple scattering CBS shows small interference con- 
trast, because ladder and CBS contributions probe dif- 
ferent field correlations, and thus produce a small total 
enhancement factor. For parallel polarization channels, 
1 1| I and h |j h, the enhancement factor tends to a larger 
value. In the h\\h channel, the coherent and incoher- 
ent contributions are asymptotically equal while single 
scattering tends to vanish. For infinite detuning, we re- 
cover the predictions for isotropic dipole scatterers due 
to Ozrin 0: 2.0 (h\\h), 1.25 (h±h), 1.75 1.12 
(l±l). 

Note, however, that this limit is reached very slowly: 
the enhancement hardly exceeds 1.5 in the h\\h channel 
for the largest detuning in Fig. [21 although it is twice 
as large as the total splitting A of the hyperfine multi- 
plet. The reason for the slow recovery is simple to under- 
stand: for a semi-infinite medium, long scattering paths 
contribute significantly to the CBS cone, with all paths 
beyond order N giving an integrated contribution scaling 
like N~ x / 2 . At large detuning, the Wi coefficients tend to 
their limiting values like 5/ A. This in turns implies that 
a fraction ~ 5/ A of the perfect contrast is lost at each 
scattering event, putting an effective cutoff ~ A/<5 on 
the scattering orders contributing to the CBS cone. Al- 
together, this implies that the asymptotic value 2 for the 
enhancement factor is reached only like a/ (5/ A. In other 
words, even a small unresolved hyperfine structure may 
significantly reduce the interference contrast. This is a 
clear illustration that CBS — and any quantum interfer- 
ence for that matter — is very sensitive to small couplings 
to uncontrolled degrees of freedom. 



angular momentum parameters defining the internal res- 
onances. We have examined under which conditions the 
various resonances conspire to yield the classical model 
of an isotropic dipole. As an application of the theory to 
multiple elastic scattering from hyperfine multiplets, we 
have calculated the CBS enhancement for a hypothetic 
half-space with a homogeneous average density of atoms. 

The present paper makes two major restrictive as- 
sumptions that should be lifted in subsequent investi- 
gations. First, the results presented in Sec. IIVI for the 
CBS enhancement factors have been calculated for the 
case of a closed transition with a unique ground state as- 
suring purely elastic scattering. The more general case of 
inelastic scattering from open transitions can be treated 
along the same lines by using the analytical theory de- 
veloped in section^ Second, real experiments with cold 
atoms are not performed on semi-infinite homogeneous 
media, but on inhomogeneous clouds of finite optical 
thickness. Their mean free path and thus the optical 
thickness change drastically under detuning from reso- 
nance. This in turn modifies the relative weight of differ- 
ent scattering orders and therefore also the enhancement 
factor. Consequently, for scattering media of finite ex- 
tent, these purely geometrical effects must be taken into 
account. An efficient numerical Monte Carlo method for 
this situation has been presented in [30l | . The above ana- 
lytical expressions for the mean free path and the inten- 
sity vertex coefficients can be introduced into this pro- 
gram and permit to obtain accurate CBS enhancement 
factors and full peak shapes for quantitative comparison 
with experiments. However, both questions are beyond 
the scope of the present contribution and will be dis- 
cussed elsewhere. 



V. CONCLUSION 

In the present paper, we have considered multiple scat- 
tering of photons by a disordered medium of atom at rest 
with several internal resonances. We have developed the 
analytical calculation of the total photon scattering cross 
section ct(cj) and the differential cross section (dcr/dfi) 
as a function of the light frequency and of the relevant 
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